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Motions of fluctuating Brownian particles in an incompressible viscous fluid have been 
studied by coupled simulations of Brownian particles and host fluid. We calculated the veloc- 
ity autocorrelation functions of Brownian particles and compared them with the theoretical 
results. Extensive discussions have been made on the time scales for which our numerical 
model is valid. 

§1. Introduction 

The motions of small fluctuating particles in viscous fluids have been studied for 
a long time. Although theoretical or numerical analysis based on the coupled motions 
of the particles and the host fluid are very complicated, it becomes rather simple if 
one considers only the particles' motions by assuming that the host fluid degree of 
freedom can be safely projected out from the entire degree of freedom of the disper- 
sions. One of such models is the well-known generalized Langevin equation(GLE) 
for Brownian particles, i.e., 



f rfs^/;,(t-s)V5(s) + G ( (t), (14) 

J — oo 



{G i (t)-G j (0)) = 3k B TT lJ (t), (1-2) 

where Mj and V{ denotes the mass and the translational velocity of the i-th particle, 
respectively. Tij{t) is a friction tensor, which represents the effect of hydrodynamic 
interactions (HI) between z-th and j-th particles. Gi is the random force acting on 
the i-th particle induced by thermal fluctuations of the solvent, ks is Boltzmann 
constant, and T is the temperature of Brownian particles. 

For a single spherical particle (i = 1) immersed in a infinitely large host fluid, 
the analytic form of the time-dependent frictiorP is known as 

J dtTii{t) exp(— iut) = Ai(— iw) = &itna{l + a^f—iuj/v — iuoa 2 /9v) (T3) 

where Tn(—iuj) is the Fourier transform of -Tii(^) arid w is the angular frequency. 
The first term corresponds to the normal Stokes friction for a spherical particle of 
radius a in a Newtonian fluid whose viscosity is n. The second term represents the 
memory effect, which is related to the momentum diffusion in a viscous medium. 
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Here the kinematic viscosity is denned as v = rj/pf with pf being the density of 
the fluid. The third term corresponds to the effect of the acceleration of the host 
fluid surrounding the tagged particle when the particle is accelerated through the 
host fluid. Using Eq. (|l-3p . the hydrodynamic GLE can be solved analytically. The 
translational velocity autocorrelation function (VACF) (Vi(t)-Vi(0)) /3 then obtained 
is known to exhibit the characteristic power-law relaxation for long-time region, 
which is widely known as the "hydrodynamic long-time tail" EJHSJU'EJ 

For dispersions composed of many particles interacting via HI, the situation is 
still not straightforward because we do not know the true analytic expression for the 
hydrodynamic friction tensor r^it). Some approximated expressions, such as Oseen 
or Rotne-Prager-Yamakawa(RPY) tensor, can be obtained by introducing the Stokes 
approximation, however, those expressions completely neglect the memory effect that 
corresponds to the second term of Eq. (jl-3p . This means that the hydrodynamic long- 
time tail can not be reproduced correctly with Oseen or RPY tensor. 

In the present study, we developed a numerical method to take into account 
the effects of hydrodynamics directly by simultaneously solving the Navier-Stokes 
equation for the host fluid with the Brownian motions of the particles. We first 
examined the VACF for a single Brownian particle and compared it with the analyt- 
ical form mentioned above. Secondly, we examined the rotational motions of a single 
Brownian particle. We furthermore examined the motions of Brownian particles in 
harmonic potentials to check the validity of our method. 



§2. Simulation method 



Here we briefly explain the basic equations of our numerical model since those 
are explained in detail elsewhereP 1 A smooth profile function < (f)(x, t) < 1 
is introduced to define fluid ((/> = 0) and particle (0 = 1) domains on a regular 
Cartesian grid. Those two domains are separated by thin interface regions whose 
thickness is £. The position of the z-th particle is Ri, the translational velocity is Vi, 
and the rotational velocity is J?j. The motion of i-th particle with mass Mj and the 
moment of inertia Ij is governed by the following Langevin-type equations, 

M i d ^ = Ff + Ff + Fr + Gf, ^ = V l , (2-1) 
Ii-^ = N? + G? , (2-2) 

where F^ and Nf 1 are the hydrodynamic forces and torques acting on the i-th 
particle due to HI, respectively. Ff and Ff x denote the direct particle-particle 
interaction and external force. G\ and Gf are the random force and torque due to 
thermal fluctuations defined stochastically as 

(GY) = (Gf ) = 0, (2-3) 
(Gj(t) ■ Gj(0)) = 3k B Ta v 5(t)5 l3 , (2-4) 
<Gf (i) • Gf (0)) = 3k B Ta n 5(t)5 lv (2-5) 
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where a and a are parameters to control the temperature T. 

The motions of the host fluid are governed by the Navier-Stokes equation 

p f {d t v + v-Vv) = -Vp + T]V 2 v + pf(j)f p (2-6) 

with the incompressible condition V • v = 0, where v and p are the velocity and the 
pressure fields of the host fluid, respectively, and (f>f p is the body force defined so 
that the rigidity of the particles is automatically satisfied. Note that Ff 1 and N- 1 
are determined from the body force </>/pP® 

§3. Results and discussion 

A single spherical particle fluctuating in a Newtonian fluid was simulated in the 
absence of external forces Ff x = as depicted in Fig. 1 . We take the mesh size A and 
r = A 2 pf /rj as the units of space and time. Simulations have been performed with 
7] = 1, a = 5, and £ = 2 in a three-dimensional cubic box composed of 64 x 64 x 64 grid 
points. The particle and fluid densities are identically set to be unity, p p = pf = 1. 

Figure 2 shows our simulation results (A) of VACF for a single Brownian par- 
ticle fluctuating in a Newtonian host fluid at ksT = 0.83. The temperature T was 
determined by comparing the long-time diffusion coefficient D s i m obtained from sim- 
ulations with D v = ksT /QT\"qaK{<P), where K(<&) takes into account the effects of 
finite volume fraction® and ^ denotes the volume fraction. The volume fraction of a 
single particle is ^ = 0.002. One finds that the VACF approaches asymptotically to 
the power-law line with the exponent —3/2, and the long-time behavior of our sim- 
ulation agrees well with the analytical solution® of the hydrodynamic GLE rather 
than the Markovian VACF which neglects memory effects. This behavior indicates 
that the memory effects are accurately taken into account. Similar to the transla- 




Fig. 1. A snapshot of a single Brownian particle immersed in a Newtonian fluid. The one eighth 
of the entire system is graphically displayed. The color map on the horizontal plane shows the 
value of the local fluid velocity in the x direction. 
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Fig. 2. The translational velocity autocorrelation function (Vi(t) • Vi(0)}/3 (triangle) and the rota- 
tional velocity autocorrelation function (f2i(t) ■ J2j(0))/3 (circle) for a single Brownian particle 
fluctuating in a Newtonian fluid. The simulation data was taken at fcsT = 0.83. The solid 
lines indicate the analytic results for the translationaP' and the rotational motions. The dot- 
ted lines show power-laws, Bt~ 3 ^ 2 with B = fcsT/12p/(7ri/) 3 for the translational motions and 
Ct~ 5 ^ 2 with C = 7rfcsr/32p/(7r!/) 5//2 for the rotational motions. The dashed lines indicate 
the Markovian VACF and RVACF, which decay exponentially as exp(— t/rs) and exp(— t/r r ), 
respectively. 



tional motions, we have studied the rotational motions of the Brownian particle in 
the host fluid. The GLE of the rotational motions for a single spherical particle can 
be written as 



hfli = - I ds^{t - s)f2i(s) + Gi(t), 

J — oo 

(Gt) = 0, (Gi(t) • Gi(0)> = 3k B Tn(t), 



(3-1) 
(3-2) 



where the time-dependent friction fi(t) has the form fi(-iuj) = 8nr]a s [l — iuo/2>v{l + 
a-\J '—iu I 'z/)p3 in Fourier space. The first term in ft is the Stokes friction and the 
second term represents the memory effect due to the kinematic viscosity of the fluid. 
The GLE can be solved analytically, and the analytical solution of the rotational 
velocity autocorrelation function (RVACF) is obtained in the following form 



(f2i(t) ■ «i(0)) 

JtksTv 
8-Krja 5 



-^exp(-yt/^) 

37T 



.3/2 



[i-(^ + l)2/] 2 + y(i-^) 2 



, (3-3) 
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where r v = a 2 jv and r r = Ij/87rrya 3 . In Fig. 2, simulation results (O) °f RVACF 
are also plotted. The RVACF clearly shows the asymptotic approach to the hy- 
drodynamic long-time tail with the exponent —5/2 which agrees well with the an- 
alytical solution (|3-3p rather than a simple Markovian RVACF. By comparing the 
present simulation results with the corresponding analytical solutions more in de- 
tail, one may notice that some discrepancies become notable for t < tb or t < r r , 
where tb = Mi/6irna = 2a 2 p p /9rj ~ 5 is the Brownian relaxation time and r r = 
Ii/8iri]a 3 = 3t b /W ~ 1.5 is the Br ownian rotational relaxation time. For opposite 
cases t > tb or t > r r , however, the agreements between the numerical results and 
the analytical solutions are excellent. This is because we neglected memory effects 
in thermal noises Gf and G^. We however believe that the long-time behavior of 
our numerical model is valid for t > t b since tb is much longer than the memory 
times of the thermal noises. 

There exist many other characteristic time-scales in particle dispersions. Impor- 
tant ones are the kinematic time-scale t v = a 2 pf/rj = 25 which measures the momen- 
tum diffusion over the particle size and the diffusion time-scale td = a 2 /D ~ 3 x 10 3 
which measures the particle diffusion over the particle size. As one can see in Fig. 2, 
the present model works quite well for the time-scales comparable to t v and td, 
while it becomes inaccurate for t < tb- 

In order to test the validity of our method for the long-time behavior of Brow- 
nian particles, we next applied the present model to simulate Brownian particles 
fluctuating in external harmonic potentials. The potentials are introduced with the 
form 

Ff x = -k(Ri - H?) = -kARi, (3-4) 

where R^ q is the ith particle's equilibrium position and k is the spring constant. 

Figure 3 shows the positional autocorrelation function (ARi(t)-ARi(O)) /3 of two 
Brownian particles in harmonic potentials whose minimum positions are separated 
by a fixed distance of 5a. The pair of particles are interacting only hydro dynamically, 
and there exists no direct interactions between them. The spring constant is set to 
k = 10, and the temperature is fc^T ~ 0.0066, which was determined by the average 
potential energy k B T = k{AR 2 )/3. The simulation results (O) 

agree well with 

the hydrodynamic analytical solution^ 1 in harmonic potentials which account for 
the effects of finite volume fraction. The analytical solution was derived by solving 
the GLE of a single Brownian particle in a harmonic potential which includes the 
modified Stokes friction £ = dirr/oK (#) . The correlation functions decay much slower 
than the Markovian relaxation functions. We also confirmed that the validity of our 
method is excellent for t > tv,. 



§4. Conclusion 

We proposed a numerical model to simulate Brownian particles fluctuating in 
Newtonian host fluids. To test the validity of the model, the translational velocity 
autocorrelation function (VACF), the rotational velocity autocorrelation function 




Fig. 3. The positional autocorrelation function (ARi(t) ■ ARj(O)} /3 (circle) for a system composed 
of two Brownian particles in harmonic potentials —k(Ri — R^ q ). The distance between their 
equilibrium positions is \R1 9 — R^l = 5a. The solid line shows the analytical solution.^ The 
dotted line shows the Markovian functions, i.e. the solution of MiVi = —6irnaVi — kRi + G;. 



(RVACF), and the positional autocorrelation function of fluctuating Brownian par- 
ticles were calculated in some simple situations for which analytical solutions were 
obtained. We compared our numerical results with the analytical solutions and found 
excellent agreements between them specially for long-time regions t > tb while some 
discrepancies were found for short time regions t < tb- This is because our model is 
designed to simulate correct long-time behaviors of Brownian particles in host fluids. 
Applications of the present method for more complicated situations are in progress. 
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